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ABSTRACT 

We have used a combination of high resolution [HST ACS-HRC, ACS-WFC, 
and WFPC2) and wide-field {ESO-WFl) observations of the galactic globular 
cluster NGC 6388 to derive its center of gravity, projected density profile, and 
central surface brightness profile. While the overall projected profiles are well 
fit by a King model with intermediate concentration (c = 1.8) and sizable core 
radius {vc = 7'.'2), a significant power law (with slope a = —0.2) deviation from 
a fiat core behavior has been detected within the inner 1". These properties 
suggest the presence of a central intermediate mass black hole. The observed 
profiles are well reproduced by a multi-mass isotropic, spherical model including 

0- 



a black hole with a mass of ~ 5.7 x 10^ Mf^. 



Subject headings: Globular clusters: individual (NGC6388); stars: evolution - 
black hole physics 



^Based on observations with the NASA/ESA HST, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5-26555. Also based on WFI observations 
collected at the European Southern Observatory, La Silla, Chile. 
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INTRODUCTION 



The surface brightness (SB) and the projected density profiles of the vast majority of 
globular clusters (GCs) are well reproduced by a family of simple models characterized by 
an extended isothermal core and a tidally trucate d envelope — the so- c alled King models 



(King 1966). However a fraction (~ 15-20%, see iDjorgovski &: King I (119861 )) of galactic 
GCs deviate significantly from this behavior. The projected density profiles of these clusters 
do not exhibit an extended core, showing instead a power law behavior S(r) oc with a 
ranging from —0.8 to —1.0 . This feature has been thought to arise from the dynamical 
evolution of stellar systems that have experienced the collapse of the core. These are called 
post-core collapse clusters, hereafter PCC). 

However, other processes can affect the shape of a star cluster density profile; among 
these the existence of an Intermediate Mass Black Hole (IMBH) in the central region has 
recently received attention. Interestingly enough, detailed collisional A^-body simulations 
(Baumgardt et al., 2005, hereafter BMH05; Trenti et al. 2007) and theoretical arguments 
( iHeggie et al.ll2007l ) have shown that the presence of a IMBH yields quite a different SB profile 
than core collapse does. These studies indicate that in an initially dense cluster a IMBH 
gives rise to a strong expansion of the central region that, in turn, leads to a quasi-steady 
configuration resembling that of a medium concentration cluster with a core-like profile. 
Thus, the clusters most likely to harbor IMBHs are those having the appearance of normal 
King model clusters except in the very central regions where a power law deviation from a 
fiat behavior is expected. The exponent of this power law is predicted to be significantly 
lower [a ~ —0.2) than in the PCC case (BMH05, Miocchi 2007). Small departures from 
a King model have been observed in the projected density profile of a few GCs (included 
NGC6388) by Noyola & Gebhardt (2006, hereafter NG06). 

The confirmation of the existence of IMBHs and an estimate of their frequency would 
be important for a number of astrophysical problems like the formation processes of super 
massive BH in galax ies, super-Eddington X-ray sources in extragalactic globu l ar clusters 



(ISivakoff et al.l l2007l ) , the origin of ultraluminous X-ray sources ( iMillerl l2003l : iFabbiano 
2006h . Despite of their potential importance, the existence of IMBH in GCs is still an open 



questi on. For instance, the evid ence for an IMBH in M15 reported by Ivan der Marel et al 



see 



( 20021) and lGerssen et all (120021 ) has been questioned bv lBaumgardt et al I (l2003af). iBaumgardt et al 
(l2003bh also question evidence fo r an IMBH in Gl in M31 JCebhardt. Rich fc Ho Il2002h (but 



Gebhardt. Rich fc Ho I (120051 ) and the recent findings by Ulvestad et al (2007) and Green 



& Ho, 2007). 

Here we present accurate surface density and SB profiles obtained with a combination of 
high-resolution and wide-field observations of the galactic globular cluster NGC 6388, which 
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a number of authors (BMH05; NG06; Drukier & Bailyn 2003, hereafter DB03; Miocchi 2007) 
have indicated as a prime candidate to harbor an IMBH. These profiles nicely match the 
theoretical models of a cluster with a 5.7 x IO^Mq BH: having an extended core, intermediate 
concentration, but also significant deviations from a flat core distribution in the innermost 
cluster regions (r ^ 1"). 



2. THE DATA 



In this paper we make use of a combination of high-resolution and wide-field photo- 
metric data sets, obtained with WFPC2 and ACS on board HST, and with WFI at ESO- 
2.2 m, respectively. A detailed description of the observations, photometric reduction and 
astrometry procedures of the data obtained with WFPC2, ACS-WFC, and WFI is given in 
a companion paper discu ssing the Blue Straggler Star and Horizontal Branch populations 



(jPalessandro et al. 1120071 ). Here we use only the optical {B, V, I) samples from the entire 
multi- wavelength data set. These have been all homogenized and transformed to the John- 
son magnitude system. All the star positions have been placed on the absolute astrometric 
system using several hundred astrometric reference star s from the new astrorn etric 2- MASS 
catalogu^, following the procedure described, e.g., in Lanzoni et al. (2007), with a final 
astrometric accuracy of the order of ~ 0'.'3 both in RA and Dec. 

Additional images obtained with the ACS-HRC have been retrieved from the ESO/ST- 
ECF Science Archive. These data sample the cluster central region with a field of view 
(FoV) of 26" X 29" and a spatial resolution of 0'.'027pix~^. The HRC data were obtained 
through filters F555W (V) and FSliW (/), with total exposure times of 6 20 and 3070 s, 
respe ctively. After corrections for geometric distortions and effe ctive flux (ISirianni et al. 
2005 ). the photometric analysis was performed by using SExtractor ( Bertin &: Arnout^ 19961 ). 
adopting a fixed aperture radius of 4 pixels (O'.'IOS). The sample has then been astrometrized 
and photometrically calibrated by using the stars in common with the ACS-WFC catalog. 
The color-magnitude diagrams based on the data from all four data sets are shown in Fig. [TJ 



3. CENTER OF GRAVITY 

Given the absolute positions of individual stars in each catalog, the center of gravity, 
Cgrav, of NGC 6388 has been determined by averaging the coordinates a and 6 of all stars 



^Available at 



http : //irsa . ipac . caltech . edu 
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detected in the highest resolution data set (the HRC sample). In order to correct for spurious 
effects due to incompleteness in the very inner regions of the cluster, we considered only stars 
brighter than V = 20 (roughly corresponding to the sub-giant b r anch of the cluster). By 



following the iterative procedure described in iMontegriffo et al. I (119951 . see also Ferraro et 
al. 2004), the center of gravity is located at ajaooo = 17^ 36"^ 1?!23, (5j2ooo = -44° 44' 7'.'!, 
with an uncertainty of O'.'S in both a and 6. A careful examination of field inside the core 
radius shows that our determination of the center is biased neither by the presence of very 
bright stars, nor of a star clump. The derived Cgrav is l ocate d ~ 2'.'6 southeast (Aa = 3'.'4, 
A6 = — I'.'l) of that derived by iDjorgovski &: MeylanI (119931 ) using the surface brightness 
distribution. An accurate comparison with the center adopted by NG06 is not possible since 
the value listed in their Table 1 is just referred to the world coordinate system of a specific 
WFPC2 image, however a visual inspection suggests that their center is located ~ O'.'S SE 
of ours. 



4. PROJECTED DENSITY AND SURFACE BRIGHTNESS PROFILES 

We have determined the projected density profile of NGC 6388 using direct star counts 
over the entire cluster radial extent, from Cg^av out to ~ 1400" ~ 23'. This distance is 
significantly larger than the expected cluster extension (r^ = 372", Harris 1996). In order 
to avoid spurious effects due to possible incompleteness, only stars brighter than V = 20 
have been considered (see Dalessandro et al 2007). There are more than 58,000 stars in the 
entire (i.e., the combin ation of ACS, WFPC2 and WFI) photometric data set. Following the 



procedure described in iFerraro et al.l (jl999al ). we have divided the sample in 40 concentric 
annuli, each centered on Cgrav Each annulus has been split into an adequate number of sub- 
sectors. The number of stars lying within each sub-sector was counted, and the star density 
was obtained by dividing these values by the corresponding sub-sector areas. The stellar 
density in each annulus was then obtained as the average of the sub-sector densities, and the 
uncertainty in the average values for each annulus was estimated from the variance among 
the sub-sectors. The radius associated with each annulus is the mid-point of the radial bin. 
The outermost (r ^ 7') measures have an almost constant value and their average has been 
used to estimate the Galaxy contamination to be ~ 56 stars arcmin"^. Subtracting this 
background yields the the profiles shown in Figs. [2j 

If the innermost (r < 1") points are excluded, the density profile is well fit all over the 
entire extension by an isotropic, single-mass King model with a core of Vc = 7'.'2 and an 
intermediate concentration (c = 1.8). These values are similar to those quoted by Trager et 
al. (1995; = 7'.'4 and c = 1.7), Harris (1996; = 7'!2 and c = 1.7), and McLaughling & 
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van der Marel (2005; Vc = 7'!8 and c = 1.71)o In the inner ~ 1" the observed profile shows 
an indication of a deviation from a flat core behavior. With only 7 stars in the innermost 
bin, the statistical error of the inner bin (0" < r < 0'.'3) is relatively large. With star counts 
alone the exact amount of the deviation from the fiat core cannot be reliably estimated. 

Exploiting the exceptional high resolution of ACS-HRC images we have computed the 
SB profile by direct aperture photometry to more accurately determine the inner shape 
of the cluster profile. In the innermost region (r < 1") we have used two sets of annuli 
stepped by 0'.'3 and O'.'S, respectively. The SB values have been computed as the sum of 
the counts in each pixel, divided by the number of pixels in any given annulus. The counts 
have then been converted to a magnitude scale and calibrated by using a relation that 
links the "instrumental" magnitude to the calibrated one (obtained by performing aperture 
photometry for a number of high S/N isolated stars). The resulting SB profile for the 
innermost 10" from the center is shown in Figure [3l A steepening of the profile at r ^ 1" is 
clearly apparent, in agreement with what we found above for the surface density distribution. 
A linear fit to the inner points suggest that the slope of the profile is 0.6 ± 0.06 in the 
logE — logr plane. In terms of the surface luminosity density /(r), if / oc r° we find 
a ~ —0.23 . This is steeper (but still marginally consistent within the errors) than the slope 
a = —0.13 ± 0.07 derived from the analysis of WFPC2 images by NG06. The N06 profile 
is shown for comparison in Figure 3, as can be seen their profile is fully compatible with 
our data in the common region. The use of high resolution images allow us to probe the 
innermost region of the cluster where most of the deviation from a fiat behavior occur^. 



5. DISCUSSION 

The properties of the projected density and SB profiles derived in the previous section 
for NGC 6388 are not those of a cluster which has experienced the core collapse. Instead they 
are just what BMH05 suggest as the signatures of a cluster harboring an IMBH in its center: 
(i) a. typical King profile with intermediate concentration (c = 1.8) in the external regions, 
(a) a sizable core, and (in) a inner logarithmic slope a ~ —0.2. These features have been 
recently confirmed by the predictions of a self-consistent parametric model that includes the 



■^Note that the higher concentration quoted in the present paper imphes a 20%) larger tidal radius 
(rt — 454") than previously determined. 

"^Even the small difference in the center determination can play a role. Simulations have shown that 
even an offset of only a few ff.'l is sufficient to flatten the profile. An additional difference in the slope 
determination might arise from the different approach used by NG06, who removed the bright stars and did 
not fit a power law to the data but instead calculated the derivative of the smooth central profile. 
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presence of a central IMBH ( iMiocchi 1120071 ). This model consists of a multi-mass isotropic, 
spherical King model, which has been extended inside t he region of gravitational influence 
of the BH via the inclusion of the iBahcall fc WoliI (119761 ) phase-space distribution function. 



In order to further support the case for an IMBH in the center of NGC 6388, we 
have used this model to reproduce the observed density and SB profiles. A Salpeter mass 
function (dA^ oc m~^ '^^dlogm) is assumed in the model and seven discrete mass bins are 
used to approximate the continuum mass spectrum of the real cluster. The stellar mass 
range in the interval from 0.3 to 0.9 Mq (here and thereafter, where not specified, masses 
are measured in Mq) equally subdivided in 6 bins 0.1 wide. These were populated with 
main sequence stars. In addition the bins [0.5, 0.6] and [0.7, 0.8] include WD populations 
with mass 0.55 and 0.75, coming from progenitors with mass, respectively, in [0.9, 1.5] and 
[1.5,4] ranges. The seventh mass bin contains a massive WD population with m = 1.2, 
hypothesized as remnants of stars with mass 4-8. The [0.8, 0.9] bin has (m) = 0.84 and 
contains the TO stars, plus giants and HB stars. The light-to-mass ratios were taken to 
be {4.9 X 10"^, 10~^, 2.3 x 10^^, 6.5 x 10"^, 0.19, 10, 0} corresponding to the bins ordered in 
increasing average mass. The velocity dispersion of the seven components is constrained by 
the requirernent o f complete ener gy equipartition at the border of the BH influence region 



[see iMiocchil l2006l : iMiocchi 1120071 . for details), where the adimensionalised potential Wbk, 
along with the the ratio between the BH mass (Mbh) and the cluster mass (M), determine 
the form of the generated profiles. Besides of the various scale parameters, Wbk and Mbh 
are adjusted to obtain the best fit to the observed profiles. To do this, we conservatively 
include only data from the central 100" in order to avoid possible spurious effects which 
might affect the most external points of the SB profiles because of the field contamination 
subtraction. The best fit to the observed SB profile is then found for Vc = 7'.'2, c = 1.8, 
Wbk = 11.5 and Mbh = 2.2 x lO'^M (yielding P{x^ > xlt) > 99%). The level of confidence 
remains above 97% for an IMBH mass in the range 2.1-2.4 x 10~^M. The and c values 
are consistent with the value deduced by the parametric fit of the single-mass King model 
mentioned above. The results of this parametric fit to the projected density and the SB 
profiles are shown in Figs. [2] and [3] respectively. 



By assuming the total c luster luminosity Vt = 6 .72 (IHarris 1 119961 ). and the distance 
modulus (m — M)y = 16.59 jDalessandro et al. ll2007l lR. we estimate a total mass of 2.6 x 
IO^Mq for NGC 6388, yielding Mbh ^ 5.7 x IO^MqH 



'''The distance has been obtained differentially with respect to 47 Tucanae, by assuming the distance scale 
by Ferraro et al. (1999b). 

^Note that this mass is well within the range of values (2.5-10 x 10'^) derived from the Mbh — cr relation 
(Ferrarese & Merritt 2000; Gebhardt et al. 2000; see for example equation (10) of DB03), by assuming the 



- 7- 



While the central IMBH is a possible explanation of the shape of the observed profiles, 
one might question whether this result is unique. In fact, a central concentration of mas- 
sive remnants (like white dwar fs and neutron stars) has been proposed as an alternative 
explanation in the case of M15 (Ivan den Bosch et al II2006I ). However, we have found that a 
multi-mass King model including a population of such remnants but without central BH is 
unable to reproduce the observed slope of the profiles in the core region. Since our evidence 
is based exclusively from the shape of the density profile, the presence of a IMBH at the 
center of NGC 6388 is still debatable. Accurate kinematical studies of the motion of stars 
in the central region of the cluster are needed to solidify the case. 

The region in which to seek for the possible kinematical signatures of a BH is very 
small. The self-consistent model here employed generates a projected velocity dispersion 
profile that shows a sharp rise from the "isothermal plateau" to a purely Keplerian behavior 
at r ~ &!16 (i.e., ~ 0.02 rc). A more promising path to detect the kinematic signature of a 
BH is by proper motion measurements (DB03). Some stars should move with anomalously 
high velocities under the influence of the BH. In order to estimate the number of these stars, 
we first need to evaluate the BH radius of influence r/j. A crude estimate of Vh is given by: 
Th = G MRn/cTn, where an is the velocity dispersion outside Vh. By assuming ctq = 18.9 km s^^ 

0.07 pc, corresponding to I'.'l with the cluster distance 



Prvor & Mevlan 


1993 


), we find Vh 


(d = 13.18 kpc; 


Dalessandro et al. 



20071 ). By using equation (7) of DB03, it is possible 



to estimate the number of stars (A^), measurable through proper motions studies, having 
velocities 3 times the cluster velocity dispersion outside r/^: A^ = 0.27 Sor^, where So'"^ is 
the number of stars within r^. This relation suggests that ~ 27% of stars within rh are 
expected to show anomalously high velocity. We can directly derive this number from the 
HRC images. Adopting the value of the cluster center presented above we count 28 relatively 
bright (y < 19) stars within = I'.'l (out of a total of ~ 85 stars detected down to ~ 22), 
corresponding to a total of ~ 7 high-velocity stars. This estimate shows that the presence 
of ~ 5 X IO^Mq BH in the center of NGC 6388 can in principle be kinematically confirmed 
in the near future through accurate proper motion studies or radial velocity measurements 
with Adaptive Optic supported instruments. However, these measurements can be quite 
challenging. According to Figure 1 of DB03, the high velocity stars are expected to be 
mainly confined within OAr^ (only ~ 0'.'4 from the center). They would have speeds of order 
60kms~^ which, given projection effects, would g ive tangental velocities of order 20kms~^. 
The distance derived by iDalessandro et al. I (120071 ) and adopted here is 25% larger than that 
from Harris (1996), so the resulting proper motions would be ~ 0.3masyr~^. Given the 
current estimate for the proper motion measurement error of roughly 0.3 mas, a baseline of 



low {(3 = 4.02) and the high (/? = 4.65) exponents, respectively. 
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at least 3-5 yr would be required. 

This research was supported by contract ASI-INAF 1/023/05/0, PRIN-INAF2006 and 
by the Ministero dell'Istruzione, dell'Universita e della Ricerca. ED is supported by ASI. 
RTR is partially supported by STScI grant GO- 10524. This research has made use of the 
ESO/ST-ECF Science Archive facility which is a joint collaboration of the European South- 
ern Observatory and the Space Telescope - European Coordinating Facility. 
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Fig. 1. — Color Magnitude diagrams for the four data sets used. The WFI data are domi- 
nated by field contamination. 
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Fig. 2. — Observed surface density profile obtained by star counts from the combined pho- 
tometric data-set: ACS-HRC {filled circles), WFPC2 (triangles), ACS-WFC (squares), and 
WFI (empty circles). The background level (see Sect. H]) has been subtracted. The single- 
mass King model that best fits the profile excluding the innermost (r ^ 1") points is shown 
as a dashed line. It is characterized by a sizable core radius (vc = 7'!2), and an intermediate 
concentration (c = 1.8). The solid line shows the profile of the model (see Sect. ED including 
a 5.7 X IO^Mq BH in the cluster center. The profile from Trager et al (1995) is shown (dotted 
line for comparison. 
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Fig. 3. — Surface brightness profile derived from the ACS-HRC images within 10" from the 
cluster center. Dots refer to a radial binning of C'S, triangles sample the inner 1" by steps 
of 0'.'3. The solid line shows the profile of the model obtained by including a 5.7 X lO^M© 
BH (the same as in Fig. [2]). The short- dashed line corresponds to the King model shown in 
Fig. [21 P = 0.6 is the slope of the linear best fit (see the long-dashed line) to the innermost 
points. The profile obtained by NG06 is also shown (see the dotted line for comparison. 



